Subharmonic gap structure in d-wave superconductors 
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We present a self-consistent theory of current- voltage characteristics of d-wave/d-wave contacts 
at arbitrary transparency. In particular, we address the open problem of the observation of subhar- 
monic gap structure (SGS) in cuprate junctions. Our analysis shows that: (i) the SGS is possible in 
d-wave superconductors, (ii) the existence of bound states within the gap results in an even-odd effect 
in the SGS, (iii) elastic scattering mechanisms, like impurities or surface roughness, may suppress 
the SGS, and (iv) in the presence of a magnetic field the Doppler shift of the Andreev bound states 
leads to the splitting of the SGS, which is an unambiguous fingerprint of d-wave superconductivity. 
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Although the d-wave scenario is emerging as the new 
paradigm in superconducting cupratesotffl, there is still 
a lacking consensus with regards to current transport in 
SIS junctionsu. By now it is clear that tunneling in un- 
conventional superconductors is a phase-sensitive tech- 
nique. An important consequence of an order parameter 
showing a sign change in different momentum directions 
is the formation of Andseev bound states at zero energy 
confined to the surfacetffl. One of the consequences of 
the presence of these bound states is the appearance 
of a zero-bias conductance peak in the current-voltage 
characteristics. This peak has been observed in differ- 
ent experiments using SIN-type junctions!] and SIS grain 
boundaries!! However, there is still a controversy re- 
lated to the observation of the subharmonic gap struc- 
ture (SGS) in cuprate SIS junctions. Devereaux and 
Fulde suggested that the Andreev scattering in SNS con- 
tacts could, be used to identify the symmetry of the order 
parametertl. Based on their analysis the SGS revealed 
in different experiments in YBCO contacts were inter- 
preted as an evidence of the existence of a well-defined 
nonvanishing gapE3lIll. After more detailed but non-self 
consistent analysis, different authors have concluded that 
the SGS are weak in d-wave superconductors, since av- 
eraging over .the anisotropic gap washes out any promi- 
nent featureaijO. However, there are also more recent 
experimental, reports of SGS in YBCO edge Josephson 
junctionstilE-3, and eiien in the c-axis tunneling of differ- 
ent cuprate contactsE3. Auerbach and Altman proposed 
an alternative interpretation of the appearance of pro- 
nounced SGS in Ref. claiming that this structure is 
an indication of magnon tunnclingrthat can be explained 
in the context of the SO(5) theoryO. 

In this Communication we present a self-consistent the- 
oretical analysis of the ac Josephson effect in d-wave su- 
perconductors. Using quasiclassical methods we deter- 
mine the local electronic properties of the superconduc- 
tors in the interface region. This includes effects on the 
order parameter (OP) profile and on the local density of 
states (DOS) by pair breaking caused both by quasipar- 
ticle scattering off the interface andj-afl 
distributed impurities in the crystalsll§E £ 
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I-V characteristics of the d-wave/d-wave SIS-junction are 
computed using the local surface Green's functions con- 
taining all relevant information of superconducting elec- 
trodes and solvingJie appropriate boundary conditions 
for a point contacts. We shall show that, not only does 
the SGS survive after averaging for all crystal misorien- 
tations, but there are qualitatively new features in the 
I-V characteristics that can only be captured by a truly 
self-consistent calculation. For instance, there is an even- 
odd effect in the SGS which originates from processes 
connecting the zero-energy states with Andreev states 
formed close to the gap edges. Such an even-odd effect 
in the differential resistance has been-xeported in YBCO 
edge-junctions by Nesher and Korer£3. These states are 
missed in a non-self consistent calculation where the sup- 
pression of the OP close to the interface is neglected. We 
also show that the SGS in d-wave contacts can be tuned 
with a magnetic-field. The Doppler shift of the Andreev 
bound statesc2lEj leads to a very peculiar splitting of the 
SGS, which provides a clear signature of d-wave super- 
conductivity. 

In the limit of not too low interface transparency, 
D > 0.1, the main feature of the I-V characteristics 
of conventional SIS contacts is the appearance of SGS, 
which consists of a series of maxima in the conductance 
at voltages eV = 2A/n, where A is the superconducting 
gap and n is an integer number. This peculiar structure 
is due to the occurrence of multiple Andreev reflect fops 
(MAR) and by now-well established both theoreticallycJ 
and experimentallyEZI in the case of s-wave superconduc- 
tors. In this work, our goal is to extend the analysis of the 
SGS to the case of superconducting cuprates. For this 
purpose, we consider a voltage biased contact, consist- 
ing of two d x 2_ y 2 superconductors separated by a single 
interface of arbitrary transparency. The order param- 
eter on side i, i = L, R, is rotated by aj with respect 
to the surface normal, and we denote junction type by 
the relative crystal orientations as d aL -d aR . There are 
several experimental realizations of this system, among 
which the. bicrystal grain-boundary junctions are ideal 
examples!! We carry out the calculation of the current 
following the formulation introduced by two of the au- 
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thors in Ref. ]24j| . and refer the reader to this work for 
all the technical details. In the case of a constant bias 
voltage, V, one can show that the current oscillates in 
time with all the harmonics of the Josephson frequency, 
i.e. I(t) = J2 m I m e m ^\ where 0(t) = <f> + (2eV/K)t 
is the time-dependent superconducting phase difference. 
We concentrate ourselves in the analysis of the dc cur- 
rent, denoted from now on as I. Furthermore, we as- 
sume that the interface conserves the momentum of the 
quasiclassical trajectories, which allows us to write the 
current as a sum over independent trajectory contri- 
butions: I = 1/^/2^ ^(ft) cos (ft)> where ft de- 
fines the Fermi surface position. For the angular depen- 
dence of the transmission coefficient we use the expres- 
sion £>(ft) — D cos 2 (ft)/[l — Dsin (ft,)], resulting from 
a (5-like potential. Here D is the transmission for the 
trajectory perpendicular to the interface. 
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FIG. 1: (a) Angle-averaged local DOS at the interface for 
a 45° misorientation (clean case): self-consistent and non- 
self consistent, (b) Local DOS for a 45° misorientation for 
different values of the bulk-impurity scattering rate F (Born 
scatterers) , measured in units of 2irTc , where Tc is the critical 
temperature in the clean case. A is the maximum bulk gap 
for the clean superconductor, (c) Angular resolved DOS for a 
trajectory pp = 45° for different magnetic fields. All the DOS 
are normalized to the normal state DOS, Nf. (d) The current 
density parallel to the interface, j y , at different applied fields 



for a d v /4-d_ 



r/4 



contact. The distance to the interface, x, is 



normalized with the clean coherence length, fo- 



Let us start by analyzing the case of a symmetric do-do 
junction in the clean limit. In this case, the order param- 
eter is constant up to the surface, and there are no bound 
states for any trajectory. The I-V characteristics of a 
single trajectory, /(ft), coincide with those of isotropic 
s-wave superconductors, and exhibit a pronounced SGS 
at eV — 2A(ft)/n pa]. Since the different trajectories 
see different gaps, the relevant question is whether the 
SGS survives after averaging. The answer can be seen 
in Fig. ^, where both the current and differential con- 
ductance are shown for different transmissions. One can 



clearly see a SGS at eV = 2A/n, although it is more 
rounded than in the s-wave case. Notice also that the 
SGS are better defined as minima in the conductance at 
the subharmonic voltages, i.e. as maxima in the differ- 
ential resistance. At this point a natural question arises: 
why is the maximum gap the energy scale revealed in the 
SGS? The idea is that the jump of the current at the 
opening of a new MAR scales with the gap of the corre- 
sponding trajectory. Therefore, the trajectories with the 
largest gaps dominate the contribution to the SGS. 
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FIG. 2: do-do contact in the clean case: (a) I-V char- 
acteristics for different transmissions, from bottom to top 
D = 0.1,0.2,0.4,0.6,0.8,0.95,1.0. (b) Differential conduc- 
tance normalized by the normal state conductance Gn- The 
dotted vertical lines indicate eV" = 2A/n. 



Let us now consider the case of a d T /4-rf_ w /4 junction. 
In this case, assuming specular quasiparticle scattering at 
the interface, an AndreeK bound state forms at zero en- 
ergy for every trajectoraj. This implies that the surface 
acts as a pair-breakeiO and the gap is depressed in the 
vicinity of the interface, vanishing exactly at the barrier. 
This order-parameter profile induces the appearance,-G|f 
new bound states at the gap edges for some trajectories^, 
which are also visible in angle-averaged DOS, see Fig. la. 
As we show in Fig. la these gap singularities are not 
present in the non-self consistent calculation. This fact 
has important consequences in the I-V characteristics. As 
is shown in Ref. |24f| , assuming a constant order parame- 
ter the most prominent feature of the I-Vs is the appear- 
ance of SGS at eV = A/n, instead of in eV — 2A/n as in 
conventional superconductors. Its origin can be under- 
stood as follows. Inside the gap the current is dominated 
by multiple Andreev reflections (MAR) . In this case there 
are two types of MAR processes: (a) those which connect 
the bound states with the gap edges and (b) the usual 
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ones connecting the gap edges. The first ones give rise 
to the SGS at eV — A/n, while the second could give 
rise to the series eV = 2A/n. However, for the non-self 
consistent calculation the DOS at the gap edges is rather 
small (see Fig. la), and thus, at the opening of this sec- 
ond type of processes their probability is small. In the 
self-consistent case the main consequence of the presence 
of the gap singularities is the recovery of the odd terms in 
the SGS series eV — 2A/n. The odd terms appear now 
due to the enhancement of the probability of the MARs 
connecting the gap edges. However, the bound states at 
the gap edges do not appear for every trajectory, which 
weakens the SGS due to these processes. In this sense, in 
Fig. H one can clearly see a difference between the even 
and odd maxima of the conductance. This even-odd ef- 
fect was reported in Ref . [^5) . This result show clearly the 
relevance of the self-consistency, even in the ideal clean 
case. 




FIG. 3: d n /4-d- n /4 contact in the clean case: (a) I-V charac- 
teristics for the same transmissions as in Fig. 2a. (b) Differ- 
ential conductance. 

In d-wave superconductors the order parameter is very 
sensitive to scattering from nonmagnetic impurities and 
surface roughness. In particular, it is known that these 
elastic scattering mechanisms provide an intrinsic -hroad- 
ening for the zero-energy bound states (ZEBS)ElEl. For 
the case of Born scatterers this broadening is oc VTA, 
where T = l/2r is the effective pair-breaking parameter 
locally at the surface. This is illustrated in Fig. lb for 
the case of bulk impurities. The natural question now 
is: what is the influence of the broadening of the ZEBS 
m the SGS? In Fig. | we show the differential conduc- 
tance of a c! 7I ./4-d_ T |./4 junction for different values of the 
bulk-impurity scattering rate. Notice that as the elas- 
tic scattering rate increases the SGS disappears, which 
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FIG. 4: Differential conductance for a d w / i -d_ 7T / i contact for 
different values of the bulk-impurity scattering rate F. Trans- 
mission rates as in Fig. 2a. 



can be understood as follows: the increase of density 
of states in the gap region enhances the probability of 
single-quasiparticle processes, producing the subsequent 
reduction of the probability of the Andreev processes, 
which in turn leads to the suppression of the SGS. 

Probably the clearest signature of d-wave supercon- 
ductivity in the SGS is its evolution with magnetic field. 
It is known that a magnetic field perpendicular to the ab- 
plane, H = Hz, leads to a Doppler shift in the continuum 
excitations given by v/ -p s , where the condensate momen- 
tum is p s = — (e/c)A(x)y,_with A the self-consistently 
determined vector potentials. As shown in Fig. lc, this 
means that the Andreev bound states are shifted to an 
energy which, in the limit of a large ratio A/£o, can be 
estimated to be £b(ft,) = (e/c)w/i?Asin^,, A being the 
afe-plane penetration depth. We shall use a natural field 
scale set by a screening current of order the bulk crit- 
ical current, Hq = cA/ev/X, which is of the order of a 
TeslaEHl. As can be seen in Fig. Id, the screening currents 
flow parallel to the interface and in opposite directions in 
both electrodes, which means that the trajectory resolved 
DOS of the left and right superconductors are shifted by 
2eb(^,) relative to each other. This shift modifies the 
threshold voltages of MARs starting and ending in dif- 
ferent electrodes, leading to the splitting of the peaks 
with an odd order n in the SGS. On the contrary, since 
the magnetic field produces a rigid shift of the spectrum, 
the threshold voltages of those MARs starting and end- 
ing in the same electrode are not modified. This means 
that the positions of the peaks with an even order n in 
the SGS remain unchanged. This fact gives rise to a rich 
SGS in the trajectory resolved current, which consists in 
conductance peaks at the following voltage positions: (a) 
[2A(jj,) ± 2eb(^,)]/n, with n odd, due to MARs connect- 
ing the gap edges of the left and right electrodes (the 
signs ± correspond to electron and hole processes), (b) 
[A(^,)±2e b (jj,)]/n, with n odd, due to MARs connecting 
the bound states and the gap edges of different electrodes, 
(c) 2A(j^,)/n, with n even, due to MARs connecting the 
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FIG. 5: Differential conductance of a clean d^u-d-^u junc- 
tion with D = 0.4 and different values of the magnetic field. 
The curves has been vertically displaced for clarity and dotted 
lines have been added to guide the eye. 



gap edges of the same electrode, (d) A(jj,)/n, with n 
even, due to MARs connecting the bound states and the 
gap edges of the same electrode. Finally, there is also 
structure at eV = 2eb(fy.)/n, with n odd, due to MARs 
connecting the bound states. 

After doing the angle average most of these features are 
still clearly visible. This is illustrated in Fig. || where we 
show the differential conductance of a d^/4-d_ n u junc- 
tion with transmission D = 0.4 for different values of 
the magnetic field. Starting at large voltages, the weak 
structure seen at 2A splits with applied field. Around 



eV = A there are both type (b) and type (c) processes 
leading to a maximum at eV = A, unaffected by the ap- 
plied field, as well as a field-shift of the dip just above 
A. The field dependence of the differential conductance 
is most clearly resolved at larger biases, eV > A/2, as 
the marks of the various processes begin to overlap at 
small bias. What is important to note is that it is still 
the bulk maximum gap that gives the characteristic en- 
ergy scale for the SGS. The effect of the Doppler shift on 
the SGS or the differential conductance is only prominent 
in junctions with a sizable misorientation. For junctions 
close to the do-do case, the main contribution to the SGS 
comes from trajectories close to perpendicular incidence, 
i.e. with sin^, ~ and thus having a vanishing Doppler 
shift. 

In summary, we have presented a self-consistent anal- 
ysis of the I-V characteristics of d-wave/d-wave contacts 
at high transparencies. We have shown that it is possi- 
ble to observe SGS for all crystal misorientations. We 
predict that the presence of bound states inside the gap 
gives rise to two new qualitative effects in the SGS: (a) 
an even-odd effect, and (b) a very peculiar splitting in an 
external magnetic field. These features are unique hall- 
marks of the d-wave scenario, and we hope that our anal- 
ysis will trigger off a more detailed experimental study 
of the SGS in cuprate Josephson junctions. 

This work has been supported by the EU TMR Net- 
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ported by the DFG, and by the Swedish research council. 
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